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Abstract 

We show that properties of the beam position monitor 
(BPM) system designed for the measurement of transverse 
and energy beam jitters can be described in terms of the 
usual accelerator physics concepts of emittance, energy 
spread, dispersions and betatron functions. Besides that, 
using the Courant-Snyder quadratic form as error estima- 
tor we introduce the scalar objective function which can be 
used as design criteria of a BPM measurement system with 
needed properties. 

INTRODUCTION 

The determination of variations in the beam position and 
in the beam energy using BPM readings is one of the stan- 
dard problems of accelerator physics. If the optical model 
of the beam line and BPM resolutions are known, the typ- 
ical choice is to let jitter parameters be a solution of the 
weighted linear least squares problem. For transversely un- 
coupled motion this least squares problem can be solved 
analytically, but the direct usage of the obtained solution 
as a tool for designing a BPM measurement system is not 
straightforward. A better understanding of the nature of 
the problem is needed. In this article, following the papers 
lUimO, we show that properties of the BPM measurement 
system can be described in terms of the usual accelerator 
physics concepts of emittance, energy spread, dispersions 
and betatron functions. Due to space limitation, we con- 
sider only the case of transversely uncoupled nondisper- 
sive beam motion and inclusion of the energy degree of 
freedom and multiple examples can be found in the papers 
cited above. 

STANDARD LEAST SQUARES SOLUTION 

We assume that the transverse particle motion is uncou- 
pled in linear approximation and use the variables z = 
(x, p x ) T for the description of the horizontal beam oscil- 
lations. As transverse jitter parameters in the point with 
the longitudinal position s = r (reconstruction point) we 
mean the difference Sz(r) — z(r) — z g {r) between pa- 
rameters of the instantaneous orbit and parameters of some 
predetermined "golden trajectory" z g = (x, p x ) T ■ 

Let us assume that we have n BPMs in our beam line 
placed at positions Si, . . . , s n and they deliver readings 
b c = (bi , . . . , 6°) for the current trajectory with pre- 
viously recorded observations for the golden orbit being 
b g = (6f , . . . , &^) T . Suppose that the difference between 
these readings Sb^ — b c — b g can be represented in the form 



5b^ = (x(si)-x(si),...,x(s n )-x(a n )) + ?, (1) 

where the random vector q = (cl , . . . , <r„) T has zero mean 
and positive definite covariance matrix 



K = diagfof, ct| 



(2) 



As usual, we find an estimate Sz q (r) for the difference orbit 
parameters 5z(r) in the presence of BPM reading errors ? 
by fitting the difference in the BPM data Sb q to the known 
optical model of the beam line, i.e. by solving weighted 
linear least squares problem. If the phase advance between 
at least two BPMs is not a multiple of 180°, then the result 
of this fit is unique and is given by the formula 

5Z-(r) = (M T (r)Vf 1 M(r))~ 1 M T (r) V' 1 ■ 5%. (3) 

The calculation of the covariance matrix of the reconstruc- 
tion errors is also standard and gives the following result 



V z {r) ^ V(6z,(r)) = (M T (r) Vf 1 M (r))" 
Here 



M = 



a n (r, si) a 12 {r, s x ) 



on(r, s n ) ai 2 (r, s„) 



(4) 



(5) 



and an(ri, r%), 012(^1 , ^2) are the elements of a two by 
two symplectic matrix A(r%, r 2 ) which transfers particle 
coordinates from the point with the longitudinal position 
s = ri to the point with the position s = r 2 . 

For the considered one-dimensional case the matrix in- 
version in the right hand side of the formula (|4} can be done 
analytically and the elements of the error covariance matrix 
V z (r) can be obtained in the explicit form as follows 



1 

m— 1 



aizir, s m ) 



(6) 



(V z (r)) 1<2 = (V z (r)) 2<1 



1 ( a ll( r > S m)\ / 012 (r, S rn ) 
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an(r, s m ) 



(7) 



(8) 
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where 



1 

k,m— 1 



Note that the error Twiss parameters can also be found 
as solution of the minimization problem 



an(r, s k )a\ 2 (r, s m ) - an(r, s m )a 12 (r, s k ) 



a k cr„ 



(9) 



In theory, the formulas ©-(O contain all information 
which one has to know in order to be able to design a BPM 
measurement system with needed properties or in order to 
be able to compare expected performance of two different 
measurement systems. In practice, unfortunately, the sim- 
ple knowledge of formulas ©-(O is quite far from being 
sufficient for these purposes. Let us assume, for example, 
that we want to compare resolutions of two BPM systems 
which are supposed to be used for trajectory jitter determi- 
nation and are installed in two different beam lines. For this 
purpose we need, at least, to have a criteria how to compare 
two covariance matrices and to know how to chose the re- 
construction points (own for each measurement system) in 
which such comparison has to be done. Does all that looks 
to be fairly straightforward? 

ERROR TWISS PARAMETERS AND 
COURANT-SNYDER QUADRATIC FORM 
AS ERROR ESTIMATOR 

An important step in solving the problems marked at the 
end of the previous section was made in (U 12, where dy- 
namics was introduced into this problem which in the be- 
ginning seemed to be static. When one changes the po- 
sition of the reconstruction point, the estimate of the jitter 
parameters propagates along the beam line exactly as a par- 
ticle trajectory and it becomes possible (for every fixed jit- 
ter values) to consider a virtual beam consisting of virtual 
particles obtained as a result of application of least squares 
reconstruction procedure to "all possible values" of BPM 
reading errors. The dynamics of the centroid of this beam 
coincides with the dynamics of the true difference orbit and 
the covariance matrix of the jitter reconstruction errors can 
be treated as the matrix of the second central moments of 
this virtual beam distribution and satisfies the usual trans- 
port equation 



V z (r 2 ) = A(r u r 2 )V z (r 1 )A T (r 1 , r 2 ). 



(10) 



Consequently, for the description of the propagation of the 
reconstruction errors along the beam line, one can use the 
accelerator physics notations and represent the error covari- 
ance matrix in the familiar form 

V(r)=e ( ^ W ~ a?(r) ^ (11) 

where /3 s (r) and a^(r) are the error Twiss parameters and 
e c = y/teW^r) = 1/VA (12) 
is the invariant error emittance. 



/3(r) , a(r) 



(13) 



Under the assumption that the phase advance between at 
least two BPMs is not a multiple of 180°, the solution of 
the problem (fT3l is unique, the minimum is reached at the 
error Twiss parameters and is equal to 2/e ? . 

Parametrization (fTTT i is an essential step ahead in under- 
standing of the structure of the matrix V z (r) in compari- 
son with the formulas l|6)-(|9). It introduces such impor- 
tant characteristic of BPM measurement system as error 
emittance and shows that balance between coordinate and 
momentum reconstruction errors in the point of interest is 
defined by the values of error Twiss parameters at this lo- 
cation. Nevertheless, it still does not give a single prop- 
erty to compare two different BPM systems. Fortunately, 
the beam dynamical point of view on the BPM measure- 
ment system naturally suggests us that in order to obtain 
the needed criteria we may simply use the Courant-Snyder 
quadratic form as an error estimator. 

Let Po, ao, 70 be the design Twiss parameters and 

I x (r, z) = 7o (r) x 2 + 2a (r) xp x + /3 (r)p 2 x (14) 

the corresponding Courant-Snyder quadratic form. Using 
this quadratic form we introduce the random variable 

1% = I x (r, 5z q {r)-8z{r)). (15) 

The mean value of this random variable is equal 

(II) = 2e,m p (/? s ,/3 ), (16) 

where 



m p (P q , Pq) = (/3 ? 7o - 2a s a + 7? A)) / 2. 



(17) 



The right hand side in ( fToT l does not depend on the posi- 
tion of the reconstruction point (is a number), character- 
izes the resolution of the BPM system not in some absolute 
units but in the relative units of beam sigmas and, there- 
fore, allows to compare properties of two completely dif- 
ferent BPM systems installed in two different beam lines 
and also can be used as scalar valued (not matrix valued) 
design criteria. 

TWO BPM CASE 

Let us consider two BPMs separated in the beam line by 
a transfer matrix A(s\, s 2 ) with a\ 2 7^ and assume that 
these BPMs deliver uncorrelated readings with rms resolu- 
tions o\ and 02 respectively. Often, when one works on 
optimization of the optics of two BPM system, one speaks 
about the desire to have the large beta functions at the BPM 
locations and the phase advance being equal or enough 
close to 90° . And that is completely right, if one will in- 
terpret this desire as a way to increase absolute value of 



the ai2 coefficient, because the error emittance of the two 
BPM measurement system is inversely proportional to it 



(cti er 2 ) / \a 



12| 



(18) 



Nevertheless, because due to formula (TToT i the figure of 
merit for the quality of BPM system is not the error emit- 
tance alone, but the product of the error emittance and the 
mismatch between the error and the design Twiss parame- 
ters (large mismatch can spoil the properties of the mea- 
surement system even for the case when the error emit- 
tance is small), one has to take additional care and compare 
design betatron functions with the error betatron functions 
which are given below 

= — K2I , a s (si) = — sign (ai 2 ) an, (19) 

O2 (72 
(J 2 (72 

/3?(s2j = — ax2 , a s («2) = sign (012) 022- (20) 

Let us note that though the error Twiss parameters de- 
pend on the ratio of BPM resolutions, the error phase ad- 
vance (phase advance defined by /3 ? ) is independent from 
this ratio and is always equal to an odd multiple of 90°. 

PERIODIC MEASUREMENT SYSTEM 

In this section we consider a measurement system con- 
structed from n identical cells assuming that we have one 
BPM per cell (identically positioned in all cells) and that 
the cell transfer matrix allows periodic beam transport with 
phase advance /i p being not a multiple of 180°. Addition- 
ally, we assume that all BPMs have the same rms resolution 
&bpm- in this situation the formula for the error emittance 
is rather simple and is given by the expression 



nJ3 p (si) 



(21) 



where (3p(s%) is the value of the periodic betatron function 
at the BPM locations and 



■MA.«-U-|i-^)' 



(22) 



is the mismatch between the error and the periodic beta- 
tron functions (even so we do not assume, in general, peri- 
odic betatron functions being the design betatron functions 
matched to our beam line). 

There is a rather widespread opinion that a periodic mea- 
surement system reaches an optimal performance when its 
design Twiss parameters are cell periodic and the cell phase 
advance is a multiple of 180° divided by n. Is that really 
so? To answer this question let us take the cell periodic 
Twiss parameters as design Twiss parameters and write 



(Ix) 



4a bp m 



(23) 
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Figure 1: Functions (/j, p ) shown for n — 2,3,4,5 
(magenta, red, green and blue curves respectively). The 
gray curve shows function (fi p ) . 

Looking at the formula ( 1221 one sees that the choice of fi p 
such that sm(nfi p ) = makes the error and the periodic 
Twiss parameters equal and brings the second multiplier 
in the right hand side of the formula ((23) to the minimal 
possible value. But, in general, it does not guarantee that 
the product of the two multipliers in (1231 is also minimized. 
So the answer is not or, more exactly, not necessary. 

To be more specific, let us consider a thin lens FODO 
cell of the length L as a basic unit of our periodic system 
and let us also assume that the BPM is placed in the "cen- 
ter" of the focusing lens. In this situation 



(I x ) = 2e q m p (p„ pp) 
where 



4 <m 

nL 



* n Op) = *oo (Mp) • m 2 p fa) , (25) 
^ocOp) = sin(/i p )/(l + sin(/x p /2)). (26) 

The functions ^ n for n — 2, 3, 4, 5 are plotted in figure 1 
together with their values in the points 



fc-(180°/n) 



1. 



. , n — 1, 



(27) 



shown as small circles at the corresponding curves. One 
sees that for all n the optimal performance of our measure- 
ment system (minimum of >]/„) is reached for the phase 
advance which is different from the multiples of 180° / n. 
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